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INTRODUCTION 


The magnetic properties of Type II superconductors have been studied 
extensively by many researchers [i] 1 and, presently, a relatively good under¬ 
standing exists of the sources of hysteresis and magnetic instabilities in these 
materials. However, new effects continue to be observed 2 [2] and, in addition, 
neither theoretical nor empirical models have been found which can fully 
describe or e^qolain more than limited amounts of published data [ 1] ., 3 

To contribute to a more complete understanding of the magnetic be¬ 
havior of Type II superconductors, an experiment was conducted to measure 
the way in which magnetic flux moves into and out of a sample under the 
influence of a changing applied field. Based on an extensive body of published 
information [ i] 4 and on prior experience [2], we expected to obtain two types 
of flux motion information. First, as the external field H applied to such a 
sample is varied slowly, a relatively slow flow of flux through the sample 


1. See, for example, Reference 1 and the references therein. 

2. M. S. Lubell and E. W. Urban have independently observed large mag¬ 
netic flux flow oscillations apparently related to the onset of instabilities. 

3. See, for example. Reference 1 and the references therein. 

4. Ibid. 


wall occurs. This flux flow into or out of the superconductor is understood to 
produce spatial variations in the induction or flux density B within the mate¬ 
rial. Knowledge of flux flow characteristics would then permit inferences to 
be made of the distribution, motion, and causes of such flux density gradients. 

Second, as H changes over a wide range, the induced gradients 
spontaneously collapse at more or less regular field intervals, producing very 
rapid flux motion. During these flux jumps, as they are called, redistribution 
of flux occurs throughout much or all of the sample volume. Knowledge of 
the velocity of the jumps and conditions for jump initiation would give informa¬ 
tion concerning the thermal and magnetic diffusion properties of the super¬ 
conductor, and the occurrence of macroscopic flux instabilities. 

The work reported here concerns a series of flux flow measurements 
and the conclusions drawn from them. A schematic of the sample and sample 
holder arrangement used for the measurements is shown in Figure 1 and an 
overall schematic of the measurement instrumentation is given in Figure 2. A 
hollow cylindrical sample of Nb25%Zr with outer radius a and inner radius b 
was placed coaxially in the bore and at the midplane of a superconducting 
solenoid which provided the applied field H and permitted H to be swept at 
controlled rates or held constant. Hall effect probes near the sample mid¬ 
plane at the outer wall and in the sample bore directly measured H and the 
bore field H* respectively, following a technique used first by Kim and 
coworkers [3]. 


2 













The rates of flux flow through the specimen walls, which we call 
4> = ^ , were measured by windings covering approximately three-fourths 
of the sample length. The single layer outer coil was wound directly on the 
specimen, while the multilayer inner coil was wound on a plastic insert shaft 
which also carried the inner Hall probe and three flux jump pickup coils. 

Descriptions of the high gain integrating digital voltmeters for measuring 

• • • 

<p and 0 and of other measurements are given in Appendix A. Flux flow 
voltages typically ranged from a few to a few hundred microvolts, depending 
primarily on the field sweep rate H. 

Six flux jump pickup coils provided radial and axial flux jump motion 
information. Three identical coils were mounted near the outer wall, one at 


the midplane and one at each end; the other set of three coils referred to 
earlier were similarly placed in the tube bore. The six jump signals, typ¬ 
ically a few volts maximum amplitude and a few milliseconds duration, were 
displayed simultaneously on an oscilloscope and photographed. A represent¬ 
ative set of pulses for a flux jump is shown in Figure 3. It can be seen that 
the jump started near the bottom of the sample at the outside and proceeded 
upward and inward. 
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BACKGROUND DISCUSSION 


To provide a clearer understanding of the measurements to be 
described, a brief summary of the present understanding of the general 
magnetic properties of Type II materials is in order. Although all super¬ 
conductors are diamagnetic and strongly oppose the entry of magnetic flux 
into their interiors, Type II superconductors permit limited flux penetration 
in the form of flux quanta or fluxoids. The fluxoids can be visualized [4] as 

filamentary current vortices with small, normal phase cores which carry 

h 

most of the field. The amount of flux in one fluxoid is </> 0 = “— = 

A e 

-15 

2.07 x 10 weber, where h and e are Planck’s constant and the elec¬ 
tronic charge, respectively. The field in a fluxoid decreases very rapidly 
with radial distance into the superconducting phase outside of the core, with 
the result that even at relatively high fields where fluxoids are densely packed, 
an appreciable portion of the volume of the material is superconducting and is 
available to carry nondissipative currents. 

If a sample of so-called ideal Type n material is placed in a slowly 
increasing field H, the following is understood to occur. At fields less than 
a lower critical field H , flux cannot enter the specimen because of the 

Meissner effect, a condition of perfect bulk diamagnetism. Above H , 

c 1 
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I 



fluxoids nucleate at the sample boundary and move freely into the interior. 
Mutual fluxoid repulsion leads to a regular distribution of vortices throughout 
the ideal sample, their area density being the induction B. H is related to 
B by 

if = — if- M , ( i) 

Mo 

where M is the magnetic moment due to the spontaneous diamagnetic cur¬ 
rents which flow at the specimen surface. A plot of -M versus H for this 
situation is shown in Figure 4a. 

As H increases further, more of the specimen volume is occupied by 

normal fluxoid cores and -M decreases. Above the upper critical field 

H , the sample bulk is completely normal, and the magnetization is zero, 
c z 

Since flux flow is unrestricted, the curve of Figure 4a is reversible; an ideal 

Type II superconductor is not hysteretic. When H < H < H , the super- 

Cl cz 

conductor is said to be in the mixed state. 

In a real, disordered Type II material, such as the NbZr alloy used 
in this study, fluxoid motion is strongly inhibited by a process known as 
pinning [ 1]. Pinning centers have been identified as defects, dislocations, 
voids, precipitate particles, and many other types of micro structural features 
which are nearly unavoidable in technical materials. An important result of 
pinning is the existence of magnetic hysteresis, indicated in Figure 4b for three 
hypothetical materials. Pinning also leads to rapid spatial variations in flux 
density in such materials, as shown schematically on the right side of Figure 5 
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(b) DISORDERED TYPE II 


Figure 4. Magnetiza 











for the cylindrical specimen arrangement indicated in the inset. If, for exam¬ 
ple, the applied field has increased from zero to H 1? B varies in some man¬ 
ner from jLx 0 Hj at r = a to zero at r - r*. Further increase of H can 
produce the succeeding induction distributions shown. 

Consider the Maxwell equation 


V 




( 2 ) 


For very slowly changing fields, D and D can be neglected. In our cylin¬ 
drical geometry, equation (2) then reduces to the scalar form 


dB 

dr 


~ Mo J 


(3) 


where J is a diamagnetic, azimuthal current density. Thus the local flux 

dB 

density gradient — is proportional to the local current density, and both 

vary with B. For example, when the applied field is H 1? J (B), defined 

c 1 

below, flows in the layer r t < r < a, as indicated in the left side of Figure 5. 

Two arguments demonstrate that as flux tries to overcome the pinning 

forces, it tends to move down the induction gradients. First, mutual fluxoid 

repulsion insures that motion will be toward regions of reduced flux density. 

1 dB 

Second, the Lorentz force F = JB = - — B —— existing between the local 

Ij fj,Q dr 

current density and the individual fluxoids comprising B can be seen to be in 
the direction of decreasing B . If one calls the net local pinning force at a 
point in the superconductor F^ , then one may argue intuitively that flux will 
be almost completely pinned until F just exceeds F [5]. One can define 

Ij P 

a critical pinning condition as existing when 


il 


(4) 



In this situation, maximum induction gradients exist and J(B) is called the 

critical current density J (B). Then the magnitude of the Lorentz force is 

J B and. if the functional behavior of F_ can be specified or inferred for a 
c P 

material, then the functional behavior of J (B) can be calculated, with the 

c 

aid of equation (4). 

On the basis of these last arguments, flux flow can be understood as 

the slow process of flux entry into or departure from a sample to establish 

the internal flux distributions, which are, in turn, governed by the critical 

current density function. Similarly, flux jumps result when the critical flux 

density gradients collapse. This leads to rapid redistribution of flux and to 

dissipation of the critical currents. For example, if a complete flux jump 

were to occur at H = H 2 of Figure 5, then at the conclusion of the jump, the 

field distribution would be as indicated by the dotted line, and J would have 

c 

vanished. Frequently in flux jump experiments it is found that the currents 
do not fully dissipate. Following such incomplete jumps, local regions of 
inhomogeneous flux distribution exist, together with the circulating currents 
which maintain them. 
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FLUX FLOW MEASUREMENTS 


Consider a flux flow measurement which starts from a condition of 

constant induction B. throughout the sample wall. The applied and bore 

fields are H. and H f . , respectively, and H. = —— B. = H\ . Based on 

3 3 J Mo 3 3 

the work of Coffey [6] and of Wipf and Lubell [7], it was anticipated that 
this uniform field condition would be, in fact, well satisfied at the conclusion 

of most flux jumps because of the complete dissipation of the currents. As 

/ 

H is then increased or decreased by a small increment, a small amount of 

flux will enter or leave a thin annular region at the outer wall with some 

measured flow velocities , respectively. If H continues to change 

at the same rate, then during the next instant flux must enter or leave a 

thicker annular layer and must fill or deplete the layer to a greater average 

induction increment than in the preceding moment. Therefore, the 

magnitude of the new measured flow rate ± 0 2 will be greater than that of 
• • 

±<Pi . Qualitatively, <j> will increase smoothly in magnitude with H, if 

* • • 

the field sweep rate H is constant. The sign of (p will be that of H , 

positive for increasing field and negative for decreasing field. The quanti- 
tative behavior of <fi will depend on the exact nature of J (B) through 
equation (3) and its integrals. 
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A typical flux flow voltage variation measured for the NbZr sample 
with continuous, constant field sweep rate is shown in Figure 6. Flux flow 
voltage was recorded versus time, or equivalently versus H, since H was 
constant. Immediately after each flux jump, which appeared as a large 
voltage pulse, 5 the flux flow voltage dropped to a small value, increased 
rapidly, then rose linearly and much more slowly. There was no observable 
change in the form of the flux flow voltage, if the field sweep was stopped 
after each jump to insure that the sample was fully cooled to the helium bath 
temperature. 

The combination of the rapid flow voltage rise, followed by the slower 
increase, was not immediately understandable. A qualitative explanation of 
this behavior was proposed later and is described below. The interesting 

m 9 

linear regions, each representing a constant flux flow rate increase <p , 
were found to vary systematically with field. In Figures 7 and 8 the slopes 

of these linear intervals are plotted as a function of the applied field at the 

• • • 

midpoint of each interval for several values of H . Figure 7 shows <p for 
increasing H . The flux flow slopes have a maximum at about 6 kilogauss, 
dip through a flat minimum, then rise smoothly. For decreasing H (Fig. 8), 
no peak appears. 


5. Note that this integrated pulse is simultaneous with, but obscures the 
details of, the type of pulse data shown in Figure 3. 
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FLUX FLOW CALCULATIONS 


Following the observation of these regular flux flow rate features, 

we sought to verify them analytically. A calculation was made of the flux 

distributions and flux flows to be expected in a cylindrical specimen. It was 

assumed that each smooth flux flow sequence would start from a uniform field 

condition B = u 0 H. . as in Figure 9. It was further assumed that when H had 
J 1 

changed by some amount H (less than that required to cause the next jump), 

a 

one of the two quasistatic situations shown schematically in Figure 9 would 
apply with 

H = H. ± H , (5) 

J a 

the upper and lower signs referring to H > 0 and H < 0, respectively. 

r . is the radial distance to which the change in B(r) has penetrated; it 
min 

depends on H. and H , and on the particular J (B) . 

j a c 

In order to guarantee that the induction distribution B(r) was indeed 
quasistatic, one had to assume that the processes of thermally activated flux 
creep, described originally by Anderson [8] and Kim et al. [3, 9], could be 
neglected. Thus if the field sweep were stopped at any level, the B(r) at 


18 


M438 



INCREASING 


ACTUAL B(r) DEPENDS 
ON J (B) MODEL 


DECREASING 

FIELD 

CASE 


mB 


Figure 9. Induction distribution and calculation notation 






that instant would cease changing, and <ft should very quickly become zero. 
Restart of the sweep at the same rate should then cause cj> to return imme¬ 
diately to the same amplitude and slope which existed when the sweep was 
stopped. On the other hand, the occurrence of appreciable flux creep after 
the sweep was stopped would cause a gradual decrease to zero and would 
allow B(r) to change somewhat, giving rise to a different amplitude and 
slope upon sweep restart. During the series of experiments, it was verified 
several times that a sweep stop, followed by several seconds wait, then 
sweep restart did, in fact, give the behavior expected for a quasistatic situa¬ 
tion in which flux creep was negligible. 

To indicate the complete functional dependencies of our problem, we 
rewrote equation (3) in the form 


dB( r, t) 
dr 


- Vo J c [B(r,t)] 


( 6 ) 


It remained only to specify a functional behavior of J and to integrate 

c 

equation (6) to obtain the variation of induction with radial position in the 

specimen. The boundary conditions are B(a) = ju 0 (H. ±H ) and 

j & 


B(r 


min 



Currently eight essentially different models for the field dependence 
of critical current density of Type II superconductors have been proposed in 
the literature. They are summarized in Table i, with references to the 
works in which they appear. In each model, & is assumed to be a 
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TABLE i. CRITICAL CURRENT DENSITY MODELS 


Model 

j 

c 

References 

Bean-London (BL) 

a 

10, ii 

Kim (K) 

a 

B 0 + B 

12 

Anderson-Friedel-Silcox (AFS) 

a 

b" 

8, 13, 14 

Yasukochi (Y) 

01 

B 1 / 2 

15 

Irie (I) 

QL 

B^ 1 

16 

Fietz (F) 

/ B \ 

oi exp 1 - — 1 + y 

17 

Goedemoed (G) 

a (n H - B) 

'^0 c2 

18 

Alden-Campbell-Coffey (ACC) 

it H - B 
p 0 c2 
a 

B 1 / 2 

19, 20, 21 


temperature dependent constant. 6 For the Kim model, Bq is a constant 
whose physical meaning is not clear. Empirically it has been found that B 0 


6. For the following reason, temperature effects were not considered in the 
calculations. Experimentally, the sample was fairly completely exposed 
to the liquid helium bath, so that if the thermal conductivity of the super¬ 
conductor were large enough and the rate of heat generation due to dis¬ 
sipative flux motion were email enough, a constant sample temperature 
equal to the bath temperature 4. 2*K could be expected. Wipf and Lubell 
[7] state that for a Nb25%Zr sample approximately of the present 
specimen size, a uniform, constant temperature condition will exist 
throughout, so long as H is less than about 10 4 G/s . Since maximum 
H was 200 G/s , the isothermal assumption was justified. 
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is of the order of a few kilogauss for NbZr and NbTi alloys and Nb 3 Sn. 


For the Irie model, y is a number which should lie in the range 1 < y < 2. 


Having obtained from the integration of equation (6) an expression for 


B(r, t) or, as is the case for some of the models, an expression for r(B,t), 


the instantaneous total flux <p ( t) in the sample can be calculated. Since 


induction B is flux density or flux per unit area, the total flux in the speci¬ 


men wall is given by an integral of B over the area of the cylinder cross 


section, 


a 

0(t) = 2 7r f B(r,t) rdr 
0 

r . 

mm a 

= 27 r f ijlq H. rdr +27 r J B(r,t) rdr . (7) 

0 J r . 

mm 

a 

0(t) = 7T jUo H. r 2 . + 27T f B(r,t) rdr . (8) 

l mm J 

J r 

mm 


The first term in equation (8) represents the total flux within the cylindrical 


region 0 < r < r . due to the field H. of the previous jump. Since r . 

min j mm 

decreases with increasing H , this term is a function of time. The second 


term is the total flux in the annular region r . <r^a in which the 

mm 


induction has changed from H.. 


If the specimen is a solid cylinder and b = 0 , equation (8) is valid 


only if r (H ) > 0 . For greater H , flux reaches the axis from all 
J mm a a 


sides and a reduced flux flow rate is required to establish a given induction 


on the axis. 
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If the specimen is hollow, equation (8) is valid only if r . > b . For 

mm 

greater H , flux must flow into and across the sample wall at an increased 
a 

rate to permit the bore field to increase smoothly. 

Our primary concern in these calculations was the determination of 
• • • 

4>(t) and cj)( t) . It would, therefore, have been possible to avoid a major 
calculation step by differentiating equation (8) directly. However, since the total 
flux in the sample as a function of field could be found from the experiments 
simply by measuring the area under the flux flow rate curves, an explicit 
expression for <p{ t) was considered valuable. Consequently, except for the 
Fietz (F) model, <p( t) was determined. It is clear that 


0( t) = 


a<£ 

3H 


dH 

dt 


and that 


dcp 

8H 


H 


(9) 


0(t) = 



dH 9<j> 

dt + at 




However it is assumed that H is constant so that 
By means of equation (9), it is seen that equation 


( 10 ) 


d± 

at 

do) 


0 cb • • 

= ~~ H vanishes. 
0H 

reduces to 


<p(t) = 


cP<b 

air 



(ii) 


In our calculations we defined the following parameters: 

H. uniform field in the sample after each jump, 


H (t) 

a 


/ 1 

units are l — kG 
\ Mo 


applied field change from H 

j 
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H(t) = H.± H (t) 
J a 


total applied field 


(* kG ) 


H= ±H 


rate (and direction) of applied field sweep 


(- -) 

\ Mo s ) 


B(r,t) 

a = + at 


flux density or induction in specimen (kG) 

proportionality constant (dimensions given 
for each model) 


</>( t) 

kt) 

kv 


total flux in specimen (weber, Wb) 

rate of change of flux or flux flow rate 
(Wb/s) 

flux flow acceleration or flux flow rate slope 
(Wb/s 2 ) 


H 


x = 


H ± H 

_J_L 


H H 

c2 c2 


dimensionless applied field 


H 


V 'H 


1 


c2 


dimensionless jump field 


p= 


B 


M«H 


c2 


dimensionless induction 


a 


sample outer radius (meter) 


r . 
mm 


minimum radius at which B has changed 
(meter). 


The upper and lower signs in H and H account for the increase or 
decrease in field for positive and negative field sweeps, while the signs in 
a similarly account for the diamagnetic nature of the superconducting cur¬ 


rents. In the tables, J , r . , B(r,t) or r(B,t), <M t), <p( t), and 

c min 
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2>(t) are shown. It should be understood that r . = r . (t) through the time 

^ mm mm 

dependence of H^, and that the three flux functions depend on H. and H^, as 
well as time. 

In Appendix B are listed the results of the flux flow calculations based 

on the eight critical current density models of Table 1. In each expression 

for 4>(t) the first term is seen to be the total flux contained within the outer 

wall due to H.; for each inter jump interval this term is independent of 

J 

time. 

Although calculated flux flow rates will not be shown explicitly here, 
it may be stated that the cp results of all of the eight models predict a smooth 
rise of the flux flow rate from zero following each flux jump. There is no 
suggestion that the sharp breaks in the curves of Figure 5 should have been 
observed. This clear contradiction between the calculations and the measured 
flow rates led one to question the assumptions made in the calculation. Wipf 
suggested that the problem might stem from the acceptance of a condition of 
field uniformity in the sample following each jump. If, in fact, large 
inhomogeneities remained, then the flux would not necessarily enter smoothly 
nor assume a well-behaved distribution. 
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STOP-HEAT FLUX FLOW MEASUREMENTS 


Since it was not realistic to try to modify the calculations to account 
for some unknown flux distributions, the experiment was modified to insure 
that the constant induction assumed to exist at the start of each inter jump 
region would indeed exist. A single layer, nichrome wire heater was wound 
noninductively over the outer flux pickup coil and on the central pickup coil 
shaft. It was then possible to heat the entire specimen to above its transition 
temperature at each jump field to permit the field throughout the wall to 
become uniform and equal to the applied field H_.. 

The measurements were then made in the following manner. Immedi¬ 
ately after each jump the field sweep was stopped. The sample was then 
heated for approximately 5 seconds at a heater power of approximately 10 
watts, which was sufficient to allow the flux to reach a quiescent, redistrib¬ 
uted condition. Finally, the specimen was allowed to cool for 10 seconds to 
insure that its temperature had returned to 4.2°K, then the field sweep was 
restarted. Typical data for this stop-heat mode of operation over one inter¬ 
jump interval are shown in Figure 10. There is a smooth rise of the flow rate 
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Figure 10. Typical flux flow volta ; 



1 division = <,80 sec = 149 G 

;e data; stop-heat mode. 



from zero as predicted by the calculations. The flow rate slope 0 is 
approximately constant. 7 These results confirm the supposition made in the 
preceding section that, following each jump, the unheated specimen contains 
induction inhomogeneities. 

Figures 11, 12, and 13 show flux flow rate slope data for the stop- 

heat measurements. The average slope in weber/sec 2 of each approximately 

linear-jump flow voltage curve is plotted versus H , the applied field at 

a 

the midpoint of the interval. Figure 11 shows the discrete data for three 
sweep rates and both sweep directions. Figure 12 gives the H = 49 G/s 
points of Figure 11 on an expanded scale. Figure 13 is a semilogarithmic 
plot in which curves summarize the data of Figure 11. The slopes are 
represented by a smoothly increasing function of average applied field (the 
low field peaks of Figure 7 do not occur). 


7. At fields above l5kG (approximately), the data have a moderate downward 
curvature; in all cases the flow rate slopes are taken from linear inter¬ 
polations between the endpoints of the smooth regions. 


28 



.46 


.42 


O 

FLUX 

FLOW 

RATE 

SLOPE 





.26 


A 



186 G/i 


X 


0 


X 


0 


.22 


.18 


.14 


k 


.10 

.06 

.02 

0 


X 0 


X 


2 



4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36 38 40 42 44 46 48 


AVERAGE APPLIED FIELD (kG) 


Figure 11. Flux flow rate slopes; stop-heat mode; 

all sweep rates. 
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Figure 13. Flux flow rate slopes; stop-heat mode; 

all sweep rates. 




COMPARISON OF CALCULATIONS AND DATA 


In order to compare the 0 formulas of Appendix B for the several 
critical current density models with the experimental results, the following 

procedure was used. Actual experimentally measured fields were used in 

• • 

solving the 0 equations: H. was taken to be the applied field observed 

J 

following each heating and recooling; H was one-half of the subsequent 

a 

constant slope interjump interval. Thus H = H. ± H was the midpoint 

J a 

field of each interjump region. The constant a was chosen for each model 

so that the calculation fit the H = 118. 5 G/s curve at point A of Figure 12. 

This fit point was chosen, since it applies to an intermediate field, where 

presumably all of the models should be reasonably applicable, and it applies 

to an intermediate value of H . For point A the parameters have the 

following experimental values: <p ~ 0. 0971 Wb/s 2 ; H = 118. 5 G/s; 

H. = 20. 10 kG; H =1.53kG; H =21.63 kG. 

J a 

• • 

The results of the 0 calculations are shown in Figures 14 and 15. 
The same 118. 5 G/s data have been replotted three times in each figure to 
separate the various curves for clarity. Curves for three values of B 0 for 
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Figure 14. Flux flow rate slope calculations. 







I 


the Kim model are shown (B 0 = 0 gives the Anderson-Friedel-Silcox (AFS) 

result) . For the Irie calculation, which reduces to the Bean-London, 

Yasukochi, and AFS curves for y - 1, 3/2, and 2 respectively, only the curves 

for y = 2. 5 and 3 are shown. Note, however, that from Irie T s original 

work there is no immediately apparent physical explanation for y > 2. For 

the Goedemoed and the Alden-Campbell-Coffey (ACC) calculations, H _ was 

c2 

taken as 63.6 kG, a measured value discussed in Appendix A. 

The <p computation for the Fietz model with its three adjustable 
parameters was made as follows. It was required that J c go to zero at » 
a necessary condition which Fietz et al. did not use in their work. This 

( m 0 h c2 \ 

gave the result that y = - o> exp I - —-- ) , leaving two parameters ot 

and (3. In their paper Fietz et al. empirically found (3 to be of the order of 

5 kG for two coils wound of Nb25%Zr wire. Therefore, j3 = 5 kG was chosen 

and the computation was fit at point A as usual, giving the strange behavior of 

Figure 15. Larger values of j3 were then chosen, in order to obtain a fit to 

the data at point B of Figure 12, as well as at point A. The very good fit 

in Figure 15 resulted when (3 = 25 kG . This value is so far from those 

measured by Fietz that the agreement must be considered as physically 

inexplicable as is the agreement of the Irie computation for y > 2 . 

The apparent variation of the Bean-London (BL) curve with H is 

spurious. The BL expression for 0 depends on H , but not on H., so 

^ J 

the experimental variation of the size of the jump intervals leads to an 
apparent variation of </>. This effect undoubtedly affects the curves of the 
other models but to a much smaller degree. 
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If the Fietz model is disregarded because of its unrealistic value of 

, ACC results in the best empirical fit over the widest field range. Never- 

• • 

theless, the calculated (p is 30 percent high at the lowest field jump and 
about 8 percent low at the highest field jump. 
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NEW CRITICAL CURRENT DENSITY MODEL 


The unsatisfactory agreement with the present experiment of the 

flux flow calculations for all of the eight published critical current density 

functions (Table 1) led us to seek a better functional form for J . Since 

c 

it appeared that a probable source of the failure of the models might be the 
assumption that a is independent of field, the variation of a with B was 
determined explicitly. Choosing the Kim model as a starting point, the 
successive values of a required to match the 4> formula to the series of 
measured interjump flow voltage variations were calculated. These a's 
are plotted against average applied field in Figure 16. By means of a least 
squares curve fit they are shown to lie approximately on a straight line which 
approaches zero at about 61 kilogauss. Because this maximum field is of 
the order of magnitude of H , it appeared reasonable to assume that 
ce(B) for the Kim model varies as 


I 


co 

oo 



Figure 16. 
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of Kim model a with field. 


where a is a field independent, temperature dependent constant. Inserting 
c 

equation (12) into the Kim function of Table 1 gives 


J = a 
c 


H - B 
c2 


c Bq + B 


(13) 


a functional form which is seen to vanish at the upper critical field H and 

h 02 

to approach a large, finite value a ——- at B = 0. 

C B 0 

The proof of the validity of this model lay, of course, in its ability to 
explain the present flux flow measurements and to agree at least reasonably 
well with other published results. Therefore, the total flux, flux flow, and 
flux flow slope calculations were made from equations (6), (8), (9), and 
(13). The results are summarized in Appendix C. H was given its 

C Li 

measured value of 63. 6 kG (Appendix A). Here, we have two adjustable 


parameters, a and B 0 , which were chosen to give exact empirical agree- 

c 


ment of the <fi calculation with points A and B of Figure 11; for this fit, 

= 3. 3 x 10 3 A/cm 2 and B 0 = 6. 5 kG. <j> curves were then calculated 
for all of the flux flow slope data of Figures 11, 12, and 13. The results 
shown in Figures 17 and 18 agree remarkably well with the measurements. 

We may, therefore, conclude that equation (13) is a very good representa¬ 
tion of the field dependence of the critical current density, at least for 
Nb25%Zr. It remains, of course, to determine the validity of equation (13) 
for other Type II materials. 

It is possible to verify the choices of oi and B 0 made from the flow 

c 

slope considerations by means of an independent experimental measurement. 

In Figure 5 it is observed that when H = H , the flux distribution B (r) just 

s s 
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Figure 18. Flux flow rate slope calculation; 



reaches the tube bore. It follows that r . = b at that field. Subsequent 

mm 

increase of H causes flux to enter the bore and H ! begins to increase. 

During the flux flow experiments, a similar sequence occurs each time the 

flux front starts from tL and reaches the tube bore without a flux jump 

having previously disturbed the B distribution. Figure 19 shows an 

experimental H ! versus H curve taken for H = 20 G/s. At points E, G, 

and I, the flux front clearly reached the bore and H T increased, if only 

briefly. From the expression for r(B) in Appendix C and using the original 

values of H . ot , and B 0 , we calculated the radial variation of B for 
c2 c 

selected values of H and H to match the fields on the increasing field 

j a 

portion of Figure 19. The B results are shown in Figure 20. Note that the 

calculated B curves reach the bore at the exact fields observed in the 

experiment. This confirms the empirical choice of a and B 0 from the 

c 

fit to the flux flow slope results. Note also that flux jump number 2 initiates 
just as the flux front reaches the bore wall. This last observation directly 
confirms earlier suggestions [5] that a hollow superconductor is especially 
susceptible to flux jumps just when a sudden increase of flux flow must occur 
to permit the entire bore to begin to fill with additional field. Heating in the 
specimen, resulting from the more rapid flow, reduces the pinning forces 
and a jump may ensue. 


42 


■ ■ ■■ 


ii ii i ■ m 


■ ■■ ■ ■IIIHM | 








-*-TUBE BORE-►-«- TUBE WALL-► 


Figure 20. Induction in tube wall for experiment of Figure 19; calculated from new model. 








Thus far, solutions to the very complicated formulas for <p 

given in Appendices B and C have been obtained. It is instructive to consider 

the behavior of the critical current functions which served as the starting 

point of each flux flow calculation. Using the values of the various parameters 

which gave the curves of Figures 14 and 15, the curves for all of the 

models except that of Irie are plotted in Figures 21 and 23. If J for the 

c 

new model is taken as the correct value for the present sample, the deviations 

of the different J models can be seen to reflect the deviations from the 

c 

experimental results in Figures 14 and 15. 
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Figure 22. Critical current density functions with proper 

high field behavior. 
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DISCUSS ION OF THE NEW MODEL 

Equation (13) was found empirically and a physical basis for it remains 
to be established. In particular the significance of the constant Bp is as 
obscure for the new model as it has been since it was originally introduced 
by Kim [9] (Table 1). It is striking, though perhaps only fortuitous, that 
the value of B 0 is just that field at which the flux flow rates for the con¬ 
tinuous field sweep have a relative maximum (Fig. 7). 

It is interesting to compare the Goedemoed and Alden-Campbell-Coffey 
critical current formulas with the new one (N). At low fields where 

PqH <& « mH. the three current density models vary as pqH , 
c1 cz cz 

PoH B~ 1//2 , and (B 0 +B) 1 for G, ACC, and N respectively. Above 

cz cz 

H ,, which is of the order of a few hundred gauss for Nb25%Zr, the 
cl 

Lorentz force or pinning force functions then vary as pc 0 H B, /x 0 H B 1 ^ 2 , 

cz oz 

and Mo ^ c 2 ^ b ~ +"b ) res P ec ^ ve ^y- The low field importance of B 0 in the 
new model is apparent. 

At high fields where B 0 may be neglected, the three forms may be 

summarized as J ~ (po H 0 ~B)B , where n = 0, i/2, i for G, ACC, 

c cz 

and N respectively. 
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The n dependence can be seen in Figure 22 near H . The Lorentz force and 

Cu 

the net pinning force at high fields are then 

F L = F p ~ (m ° H c2 " B)Bl " n ; B» B 0 , (14) 

and for N they vary as - B. Coffey [21] has calculated that the 

c2 

fraction of the specimen cross section which is outside the fluxoid cores and 

remains superconducting is also proportional to ju 0 H - B. Consequently 

c u 

only the new model is seen to be consistent with the quite reasonable argument 

4 

that the net pinning force on the fluxoids varies directly with the local area 
density of superconducting phase in the material. 
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SUMMARY 


In this report we have discussed new observations of magnetic 
flux flow behavior in the Type II superconductor Nb25%Zr. The rate of 
flux flow <i> into or out of a cylindrical specimen increases linearly with 
field in the interval between flux jumps, provided the process in each interval 
starts from a condition of constant induction throughout the specimen. The 

ft 4 

slopes of the linearly increasing flow rates, which represent $, are 

themselves increasing functions of average applied field. Since the flux in 

the sample and its derivatives with respect to time can be calculated from 

an assumed variation of critical current density J^(B), it is possible to study 

the applicability of various functional forms of J which have been proposed 

c 

in the literature. Although some of the eight previously advanced J models 

c 

were able to explain the present data only over limited ranges of applied 
field, a new empirical model proposed and discussed in this report is 


shown to agree with the experimental results at all fields, except near H 

» H o 2 - B ' 

where anomalous behavior occurs. The new model, J = ol —~r-rr— 

c c B 0 + B 

is well behaved at very low fields and vanishes properly at H At 

C /~i 


intermediate fields, its variation closely resembles short sample critical 
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current densities reported in the literature. Although the physical signifi¬ 
cance of the constant B 0 remains obscure, its empirical value 6.5 kG is 
close to the values found for this parameter in the simpler, related Kim 
model. It remains to provide an explanation of the new model from considera¬ 
tions of microscopic flux pinning interactions in the superconductor. 8 


8. Portions of the work discussed here are also reported in Reference 22. 


M438 
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APPENDIX A 

EXPERIMENTAL DETAILS AND ADDITIONAL MEASUREMENTS 


Because some relatively novel experimental problems were encoun¬ 
tered in this study and because the special dc amplification system used in 
the flow measurements has very good accuracy, sensitivity, and general 
utility, a more detailed discussion of the apparatus will be given. 

The flux flow measurement system employed for this study was in 
effect a dual, rapid, high gain integrating digital voltmeter with simultaneous 
digital and analog output records. One channel was used for the outer flux 
flow pickup coil, the other for the inner coil; one channel is described. All 
of the apparatus in this dc amplifier was standard Hewlett-Packard instru¬ 
mentation, but its use in the present system combination is noteworthy. A 
2212A Voltage-to-Frequency Converter (VFC) converted the low level dc 
signal into a pulse train whose repetition frequency was directly proportional 
to the input voltage. The VFC sensitivity was 10 Hz per /xV. The VFC 
pulse frequency output was counted by a 532 IB Frequency Counter at a rate 
of 10 frequency readings per second. These frequencies corresponded to the 
average input voltage in microvolts during the counting time. They were 
digitized and recorded on one of two channels of a 5050B Digital Printer. A 
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temporary data storage feature in the printer reduced the time during which 
the printer inhibited the counter to 100 ps per print cycle. Consequently 
the system monitored the input signal 99.90 percent of the time and the 
probability was very small that any significant data fluctuations, including 
flux jumps, would be missed. 

As they were recorded digitally, the frequency data were also fed to 
a 581A Digital-to-Analog Converter (DAC) which generated a dc voltage pro¬ 
portional to the number represented by the three least significant digits of 
each counter reading. This voltage was recorded on a 7100A Strip Recorder 
whose full scale deflection was adjusted to match the DAC output voltage 
range for a 0 to 999 reading. Thus, although the counter readings were often 
greater than 999, the recorder pen always remained on scale by recycling 
as often as necessary. Note, for example, the 636 (jV and 1778 pV deflec¬ 
tions in Figure 10. It can be seen that this method effectively expanded the 
width of the final analog recording, while maintaining the same absolute 
sensitivity. Ambiguities introduced by pen recycling were removed later by 
reference to the digital record. Introduction of the 0. 1 second integration 
time into the data produced small, unobtrusive steps in the final analog re¬ 
cord, as was seen in Figures 6 and 10. 

To eliminate thermoelectric voltage problems in the flux flow 
measurements, the pickup coils were formed of a continuous conductor 
from one room temperature VFC input through the coil winding in the liquid 
helium and back to the other VFC input. The VFC-to-coil connections 
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themselves were enclosed in a cardboard box to protect them from rapid 
temperature changes. Drifts in the entire system were very slow and could 
easily be corrected for by occasionally rezeroing the VFC’s. 

Although the field sweep rates were automatically and very accurately 
controlled at intermediate magnetic fields by the magnet power supply, 
problems were encountered at low (< 8 kG) and very high (> 36 kG) fields. 
In these regions for H > 0, the magnet sweep rate changed in a slow, 
reproducible manner because of changing magnetic inductance. These 
anomalous sweep rates were independently measured by means of a bare 
copper coil in the solenoid bore and the data plotted in Figures 11 through 14 
and 16 through 17 were corrected accordingly. The data in these figures 
accurately reflect the flow slope behavior which would be observed for the 
stated constant H values. 

Since one parameter of three of the important critical current models 

considered in this report is the upper bulk critical field H , its direct 

c2 

measurement was necessary. Published values [23] of H for Nb25%Zr 

c z 

are spread around 70 kG and have large uncertainties. Since H is a 

cz 

function of composition, and the exact composition for this specimen was not 

known, an independent critical field determination was made. Under normal 

circumstances H is the field at which the sample becomes normal and the 

c2 

difference H r - H between the bore and applied fields becomes zero. Since 
both branches of the critical state curve converge to H’ = H at H , the 
upper critical field is observed as the field at which the two branches merge 
with the H r = H Hall probe characteristic. 
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The NbZr sample studied exhibits an interesting property known as 


the peak effect. It is common to many Type n materials and has been attrib¬ 
uted [24] to increasing resistance to mutual interpenetration by the flux lines 

in the sample near H ^. Its existence necessitates the definition of two 

c2 

upper critical fields. High field H' versus H data are shown in Figure A-l 

in which the peak effect is quite evident; at several fields on both critical 

state curves the slow sweep was stopped and the sample heated to above T to 

c 

define the Hall probe characteristic. 

m 

Above approximately 55 kG for H > 0 the field difference H' - H, 
instead of continuing to converge monotonically toward zero, begins to 
increase with field. It passes through a maximum and then falls rather 
rapidly to zero. For H < 0, a similar maximum in H' - H is seen. This 
rapid high field convergence of the critical state envelopes permits relatively 
precise determination of the true bulk upper critical field, which will be 
simply referred to as H One would expect that any reasonably simple 
physical model which explains the magnetization and flux flow data at fields 
below 55 kG should extrapolate smoothly into the region in which the peak 
effect is observed. It follows that the value of upper critical field which must 
be used in such a model is that field toward which the magnetization is 


converging when the peak effect commences. This effective bulk upper 


critical field will be called e ff\ • 

To allow estimation of H and H . ... the X-Y recorder pen 

c2\6it; 

separation at constant H was plotted from approximately 38 kG to H 
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Figure A-i. X-Y record of bore and applied fields 


as seen in Figure A-2. The peak effect is very pronounced and at its high 
field side drops steeply and smoothly to an intercept at 69.4 ± 0. 1 kG. A 
very linear decrease occurs prior to the onset of the peak effect rise. A 
linear least squares curve fit leads to an extrapolated effective value ^ c 2 ( e ff) 
of 63.6 ± 0.5 kG, the value which was used above in the fit of the critical 
current density models to the flux flow data. 
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APPENDIX B 

CALCULATION RESULTS - PUBLISHED MODELS 
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APPENDIX C 


CALCULATION RESULTS - NEW MODEL 
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where 


A = 6 term constant (function of x, /3 0 ) 

C = 4 term constant (function of x., /Sq) 
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